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Abstract 

In the framework of kr ~factorization approach, the production and polarization of prompt 
if(2S) mesons in pp collisions at the LHC energies is studied. Our consideration is based 
on the non-relativistic QCD formalism for bound states and off-shell amplitudes for hard 
partonic subprocesses. The transverse momentum dependent (unintegrated) gluon densities 
in a proton were derived from Ciafaloni-Catani-Fiorani-Marchesini evolution equation or, 
alternatively, were chosen in accordance with Kimber-Martin-Ryskin prescription. The non- 
perturbative color-octet matrix elements were first deduced from the fits to the latest CMS 
data on ip(2S) transverse momentum distributions and then applied to describe the ATLAS 
and LHCb data on if(2S) production and polarization at y/s = 7 TeV. We perform the 
estimation of polarization parameters Xg, X$ and A^ which determine if(2S) spin density 
matrix and demonstrate that taking into account the off-shcllness of initial gluons in the 
color-octet contributions leads to unpolarized if (2S) production at high transverse momenta, 
that is in qualitative agreement with the LHC data. 
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1 Introduction 


The production of quarkonium states in high energy hadronic collisions is under intense 
theoretical and experimental study since two decades ago, when the measurements 

of prompt J/ip and T production cross sections at the Tevatron revealed a more than one 
order-of-magnitude discrepancy with theoretical expectations obtained in the framework of 
the color singlet model [J]. This fact had induced extensive theoretical activity mainly 
connected with modeling the formation of quarkonium states from unbound heavy quark 
pairs produced in hard interation. There exist two competing theoretical approaches known 
in the literature under the names of color-singlet (CS) and color-octet (CO) models. In 
general, a quark-antiquark pair is produced in a state 2S + l L <y f > with spin S, orbital angular 
momentum L , total angular momentum J and color a, which can be either identical to the 
respective quantum numbers of the resulting quarkonium (as accepted in the CS model) 
or different from those. In the latter case, the heavy quark pair transforms into physical 
quarkonium state by means of soft (nonperturbative) gluon radiation, as considered in the 
formalism of non-relativistic Quantum Chromodynamics (NRQCD) [5)16]. The quarkonium 
formation probability is then determined by the respective nonperturbative matrix elements 
(NMEs), which are assumed to be universal (process-independent) and not depending on 
the quarkonium momentum. Though not strictly calculable within the theory, the NMEs 
are assumed to obey certain hierarchy in powers of the relative quark velocity v. To the 
leading order in v, an S'-wave vector meson such as i[)(2S) can be formed from a quark pair 
produced as color singlet 3 S'j 1 ' ) or via one of intermediate color octet states 1 Sq’\ 3 S[ 8 ^ or 
3 pj 8) with J = 0, 1 or 2. 

We know already that the CS model with leading order (LO) hard scattering matrix 
elements fails to describe the experimental data on J/-0 and T production at the Tevatron 
and LHC. Including the next-to-leading order (NLO) [7J and dominant tree-level next-to- 
next-to-leading order (NNLO*) [8] corrections to the CS mechanism significantly improves 
the description of the collider data PJ. In the NRQCD formalism, a reasonably good agree¬ 
ment with the measured quarkonia cross sections can be achieved by adjusting the NME 
values which play the role of free fitting parameters [T0HI5] . This was already demonstrated 
by comparing the calculations with the ATLAS [16], CMS d) and LHCb |T8] experimen¬ 
tal data taken at yfs = 7 TeV. However, the values of the extracted NMEs dramatically 
depend on the minimal quarkonium transverse momentum px used in the fits [19j and are 
incompatible with each other (both in size and even in sign!) when obtained from fitting the 
different sets of data. Furthermore, none of the fits is able to accommodate the polarization 
measurements [20ll2l] . The fits involving low px measurements lead to the conclusion that 
the production of S'-wave quarkonia at high px must be dominated by CO contributions with 
transverse polarization (namely, by the 3 S[ 8 ' 1 channel). The latter contradicts to the unpo¬ 
larized production seen by the CDF [221123] Collaboration at the Tevatrorjj] and CMS [2171124] 
and LHCb [21] Collaborations at the LHC. To obtain an unpolarized state it is necessary to 
either assume that quarkonium production is dominated by the scalar 1 Sq 8 ^ channel [12] or 
restrict the NRQCD fits to very high px region [19]. This problem is known as ” quarkonium 
polarization puzzle” and still far from understanding. 

lr The CDF Collaboration has measured small longitudinal polarization in prompt J/■i/j production. 
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A new solution to the polarization puzzle has been guessed in the fc^-factorization ap¬ 


proach of QCD [251125] . where studies of quarkonia production and polarization have their 


own long history (see, for example, [271 [34] and references therein). The /cy-factorization 
approach is based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [35] or Ciafaloni-Catani- 


Fiorani-Marchesini (CCFM) [36] gluon evolution equations and provides theoretical grounds 
for including the effects of initial gluon radiation and intrinsic gluon transverse 



A combination of the usual CS scheme and the kr ~factorization formalism results in reason¬ 
ably good description of the data on J/0, Xc an d T production at HERA [2711291 [301133] , 
Tevatron [281132] and LHC [34] . These results are also in good agreement with more compli¬ 


cated explicit NNLO* calculations [M] performed in the collinear approximation of QCD. 


The longitudinal polarization of directly produced quarkonia in the A/p-factorization is an 
immediate consequence of initial gluon off-shellness EH- Adding the feed-down from P -wave 
states (xc and Xb) leads to essentially unpolarized prompt J/\b and T mesons [221134] . In 


the Ay -factorizat ion approach at LO, the P -wave states are produced in 2—^1 gluon-gluon 


fusion and are expected to dominate at high transverse momenta. However, the latest LHC 
data show that the feed-down contributions Xc —> J/i^+x and Xb —> T+y do not constitute 
more than only 20 — 30% of the visible cross section at large Pt values. 

So, the production of Xc and Xb mesons requires a dedicated study which will be the 
subject of our forthcoming papers, while in the present analysis we concentrate on the 
direct mechanism and only restrict to '0(25') mesons having no contamination from higher 
states. Here we present a systematic analysis of ATLAS pH], CMS [iTi] and LHCb [18] data 
collected at yfs = 7 TeV regarding the transverse momentum distributions and polarization 
parameters Xe, X^ and Awhich describe the spin density matrix of the produced 0(2S') 


mesons. 


The outline of our paper is the following. In Section 2 we briefly recall the NRQCD 
formalism and the A/p-factorization approach. In Section 3 we perform a numerical fit to the 
latest CMS data and extract the color-octet NMEs using three different sets of transverse 
momentum dependent (TMD) gluon distributions. Later in this section we check the com¬ 
patibility of the extracted papameters with ATLAS and LHCb data on the production and 
polarization of 0(2 S) mesons. The comparison is followed by a discussion. Our conclusions 
are collected in Section 4. 

2 Theoretical framework 

We start with briefly recalling the essential calculation steps. Our consideration is based 
on the following leading-order off-shell partonic subprocesses m- 


0*(* 1 )+0*(* 2 )->0 / P^] (P)+$(*), 

g*(h) + g*(k 2 ) 0' 2 3 s[ 8 \ 3 Pf 1 (p), 


(i) 


( 2 ) 


where J — 0,1 or 2, and the four-momenta of all particles are indicated in parentheses. 


The subprocesses (1) and (2) represent the CS and CO contributions, respectively. The 


2 A detailed description and discussion of the /cT-factorization formalism can be found, for example, in 


reviews 
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corresponding production amplitudes can be obtained from the one for an unspecified cc 
state by applying the appropriate projection operators, which guarantee the proper quantum 
numbers of the cc state under consideration. These operators for the different spin and orbital 
angular momentum states read [4j: 


n 


n 

n [ 3 

3 Pj} = 


s i = 

= (Pc 


= 75 (pc + m c ) /m 1/2 , 

e(S z ) (p c + m c ) /m 1/2 , 

- m c ) e(S~) (p c + m c ) /m 3/2 , 


(3) 

(4) 

(5) 


where m = 2 m c is the mass of the considered cc state, and p c and p 5 are the four-momenta 
of the charmed quark and anti-quark. States with various projections of the spin momentum 
onto the z axis are represented by the polarization four-vector e^(S z ). 

The probability for the two quarks to form a meson depends on the bound state CS and 
hctituous CO wave functions \k^(g), where the relative four-momentum q of the quarks in 
the bound state is treated as a small quantity in the non-relativistic approximation. So, we 
represent the quark momenta as 


p c = p/2 + q, p 5 = p/2 - q. 


( 6 ) 


Then, we multiply the hard subprocess amplitude A (depending on q) by the meson wave 
function \|A a )(g) and perform integration with respect to q. The integration is done after 
expanding the amplitude A around q — 0: 


A(q)¥ a \q) = A\ q=0 ¥ a \q) + q a (8A/dq a )\ q=0 ¥ a \q) + ... 


(7) 


Since the expressions for A\ q= o and 8A/dq a \ q= o are no longer dependent on q, they may 
be factored outside the integral sign. A term-by-term integration of this series employs the 
identities [1]: 

J * [q) x/4^ [ h 




( 8 ) 

(9) 


(27t) 3 " y/An 

(2tt) 3 ^ ■** 

where VS a \x) are the radial wave functions in the coordinate representation.The first term 
in (7) contributes only to S- waves, but vanishes for P- waves. On the contrary, the second 
term contributes only to P- waves, but vanishes for S- waves. States with various projections 
of the orbital angular momentum onto the z axis are represented by the polarization four- 
vector e^(L z ). The NMEs of S'-wave states are directly related to the CS and hctituous CO 
wave functions: 


(i O* 


2S+1 t i a ) 
L J 


) = 2N C (2J+ l)|77 (a) (0)| 2 /47r, 


( 10 ) 


where N c = 3 and J = 1. A similar relation holds for if P -wave states are involved. 
The CS wave function at the origin of coordinate space is known from the measured ip(2S) 
leptonic decay width. The color-octet NMEs are extracted from experimental data and obey 
the relation 


(i O* 


3 P ( 8) 


) = (2 J + 1) (O^ 


3 p(8) 

Jl n 


), 


(ii) 
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coming from heavy quark spin symmetry at LO in v. The polarization vectors e At (S' z ) and 
e^(L z ) are defined as explicit four-vectors. In the frame where the z axis is oriented along 
the quarkonium momentum vector p M = (£7, 0 , 0 , |p|), these polarization vectors read 

e,(±l) = (0,±M,0)/V2, e„(0) = (|p|, 0, 0, E)/m. (12) 


The states with definite S z and L z are translated into states with def ini te total momentum 
J and its projection J z using the Clebsch-Gordan coefficients: 

^2, L z ]l 1 J z ) e^Sz) e u {L z ). (13) 

Sz, L z 

Further evaluation of partonic amplitudes is straightforward and is done using the algebraic 
manipulation system FORM [38]. Our results for perturbative production amplitudes squared 
and summed over polarization states agree with ones [39j. Here we only mention several 
technical points. First, according to the kr ~factorization prescription [25ll2fij . the summation 
over the incoming off-shell gluon polarizations is done using the gluon spin density matrix 
in the form e tJj e* u = k^k^/k^, where k^ is the gluon transverse momentum orthogonal to 
the beam axis. In the collinear QCD limit, when |ky| —» 0, this expression converges to the 
ordinary e^e*" = —g tJ ‘ v /2 after averaging over the azimuthal angle. Second, the ^{2S) spin 
density matrix is expressed in terms of the momenta l\ and / 2 of the decay leptons and is 
taken as 

£= 3 (lK + « - /™ 2 - (M) 

This expression is equivalent to the standard one — ~9 IX1 ’ + p 9 p v /m 2 but is better 

suited for studying the polarization observables because it gives access to the kinematic 
variables describing the orientation of the decay plane. In all other respects the evaluation 
follows the standard QCD Feynman rules. 

An important point in the NRQCD formalism is connected with the emission of soft glu¬ 
ons taking place after the hard interaction is over. It is usually assumed that the emitted soft 
gluons bring away the unwanted color and change other quantum numbers of the produced 
CO system, but do not carry any energy, thus keeping the kinematics intact. However, such 
an emission contradicts to the basic QCD property: soft gluons can never be radiated as 
they are confined. In order that the quantum numbers get changed, one needs to radiate 
a real gluon with some energy E ~ Aq CD , giving us the confidence that we do not enter 
into the confinement or perturbative domains. When considering the gluon radiation from 
3 P® and 3 P 2 8 ^ states, we rely upon the dominance of electric dipole El transitions, which 
is supported by the E835 experimental data [30]. The corresponding invariant amplitudes 
can be written as follows |41j : 


A( 3 Pl 8) ^'ijj' + g)=g 1 e ^Ds) e (C°) e (V>') e (s), 


(15) 


A( 3 PF^1/ + g)=g 2P ™<!g >) eM 


k^e 


C 9 ) 
P 


PD (g) 


(16) 


where p^°\ kj?\ and e / ) c i ; 0) are the four-momenta and polarization four- 

vectors (tensor) of corresponding particles and is the fully antisymmetric Levi-Civita 

tensor. The gluon radiation from other CO states is generated according to the phase space. 
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The cross section of 0(2S') production at high energies in the fc^-factorization approach 
is calculated as a convolution of the off-shell partonic cross sections and the TMD gluon 
densities in a proton. The contribution from the CS production mechanism can be presented 
in the following form: 

a(pp -+ 0' + X) = j k iT> V 2 )fg( x 2, k 2 ti V 2 ) x 

x | A{g* + g* -+ 0' + g)\ 2 dpldk 2 1T dk 2 2T dydy g 

where f g (x, ky, p 2 ) is the TMD gluon density, and y are the transverse momentum and 
rapidity of produced 0(25) meson, y g is the rapidity of outgoing gluon and y/s is the pp 
center-of-mass energy. The initial off-shell gluons have a fraction X\ and X 2 of the parent 
protons longitudinal momenta, non-zero transverse momenta k 1T and k 2 T (k 2 T = —k 2 T 0 0, 
k 2 T = —k 2 T 0 0) and azimuthal angles <0 and 0 2 . For the CO production we have: 


f QiTT 

a(pp 0' + X) = / - - f g (x i, k 2 T , p 2 )f g (x 2 , k 2 r , p 2 ) x 

J 1 iL/ o u x 


x\x 2 sF 

x | A(g* + g* 0')| 2 d^d^dy 


d(p i d(f> 2 

2n 2 tt 


(18) 


According to the general definition [02], the off-shell gluon flux factor in (18) is definec0 
as F — 2A 1//2 (s, kf, k%), where s — (ki + k- 2 ) 2 . The squares of the corresponding off-shell 
partonic amplitudes, as being too lengthy, are not presented there but the full C++ code is 
available on requeso The multidimensional integration have been performed by means of 
the Monte Carlo technique, using the routine VEGAS [43] . 

Numerically, we have tested several different sets of the TMD gluon densities. Two of 
them (namely, JH and AO sets) have been obtained [441145] from the CCFM equation where 
all input parameters have been fitted to describe the proton structure function F 2 (x,Q 2 ). 
Besides the CCFM-evolved gluon densities, we applied the one obtained from the Kimber- 
Martin-Ryskin (KMR) prescription [46]. The KMR approach is a formalism to construct 
the TMD quark and gluon distributions from well-known conventional ones. For the input, 
we have used leading-order Martin-Stirling-Thorn-Watt (MSTW’2008) set [47]. 


3 Numerical results 

We now are in a position to present our numerical results. First we describe our input and 
the kinematic conditions. Having the TMD gluon distributions chosen, the cross sections (17) 
and (18) depend on the renormalization and factorization scales yn and yp. We set y 2 R = 
m 2 + p? r and y 2 F = s + Q^, where Qr is the transverse momentum of the initial off-shell gluon 
pair. The choice of yn is the standard one for studying the charmonia production whereas 
the special choice of p F is connected with the CCFM evolution [441145]. Following [48], we 
set 0(2 S) mass m = 3.686 GeV, branching fraction R(0' —> p + p~) = 0.0077 and use the LO 

3 The dependence of numerical predictions on the different forms of flux factor has been studied in m- 
4 lipatov® theory, sinp.msu.ru 
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(<O* pWl^VGeV 3 

(< O* [^JVGeV 3 

(( O* [ 3 S( 8) ])/GeV 3 

(O^ [ 3 P 0 (s) ])/GeV 5 

A0 

7.04 x 10" 1 

0.0 

5.64 x 10~ 4 

3.71 x RT 3 

JH 

7.04 x RT 1 

0.0 

3.19 x 10~ 4 

7.14 x 10" 3 

KMR 

7.04 x RT 1 

8.14 x 10~ 3 

2.58 x 10" 4 

1.19 x 10" 3 

tm 

6.50 x RT 1 

7.01 x 10~ 3 

1.88 x 10" 3 

-2.08 x 10“ 3 


5.29 x RT 1 

-1.20 x 10" 4 

3.40 x 10” 3 

9.45 x 10~ 3 


Table 1: The NMEs for ip(2S) meson derived from the fit of the CMS data [T7] , The NMEs 
obtained in the NLO NRQCD fits IHJITS] are shown for comparison. 


formula for the coupling constant a s {/i 2 ) with rif = 4 quark flavours at Aqcd = 200 MeV, 
such that a s (M|) = 0.1232. 

In Table 1 we list our results for the NMEs fits obtained for three different TMD gluon 
densities. We have fitted the transverse momentum distributions for prompt -0(2 S) mesons 
measured recently by the CMS Collaboration at y's = 7 TeV ra. These measurements were 
done at moderate and high transverse momenta 10 < pt < 100 GeV, where the NRQCD 
formalism is believed to be most reliable. In contrast with [TTJ|T5], we performed the fitting 
procedure under requirement that the NME values be positive only. The color singlet NMEs 
were not fitted, but just taken from the known 0(2S') —> fi + n~ partial decay width [ 48] . 
For comparison, we also present in Table 1 two sets of NMEs, obtained within the NLO 
NRQCD in [TTHT5] . The main difference betwen Em and [I5j is in that these fits were based 
on differently selected sets of data points. 

In the fcr-factorization approach, the fitted NME values strongly depend on the choice 
of TMD gluon density. We find that the 1 Sq^ contribution is compatible with zero if the 
CCFM-evolved gluon distributions are used, but is non-negligible in the case of KMR gluons. 


For the latter, the extracted value of ( O^ 1 5 , q ) is very close to the one obtained in the 

NLO NRQCD analysis [U]. Both the NLO NRQCD fits [TT1IT5] significantly (by one order of 
magnitude) exceeds the values of (O^ 3 S'{‘ s) ), obtained with all of the TMD gluon densities, 
ft is almost consistent with estimates performed by other authors [PltoTil. In contrast with 
the results [39], our fit leads to non-zero ( O^ 


‘P 


( 8 ) 


) values. Summing up, we can conclude 
that the NME values obtained by the different authors on the basis of different data sets or 
by the same authors using different gluon densities are widely spread, that spoils the belief 
in the universality of the matrix elements. 

Now we turn to comparing our predictions with the data collected by the ATLAS [ IB] , 
CMS [I?] and LHCb [T8] Collaborations. The ATLAS Collaboration has measured prompt 
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0(2S') transverse momentum distribution at 10 < px < 100 GeV at central rapidities \y\ < 2 
[IB] . The CMS Collaboration probes the transverse momentum in the range 10 < px < 
100 GeV at \y\ < 1.2 [17] , and the LHCb Collaboration works in the kinematic range 
Pt < 16 GeV and 2 < y < 4.5 [18]. In all cases the data were obtained at y's = 7 TeV. The 
results of our calculations are shown in Figs. 1 — 3. With our sets of NMEs, we achieve 
reasonably good description of the data with any of the considered TMD gluon density. 
We observe dominance of the CO contributions in the whole px range. In particular, the 
3 S[ S) contribution dominates at the large px > 25 GeV, whereas the 1 S'o 8 ' ) channel is mostly 
important at low px values. Taken solely, the CS contributions (even incorporated with the 
fc^-factorization) are unable to describe the data. They are important at relatively low px 
only and are comparable there with the 3 Pj 8 ' ) contributions. At moderate and high transverse 
momenta the CS contributions are below the data by about one order-of-magnitude. The 
predictions obtained with the chosen TMD gluon densities are very close to each other at 
px > 6 GeV, while the difference becomes only sizable at low px < 6 GeV (see Fig. 3). 
Therefore, similar to the collinear QCD factorization, including the low px data to the fit 
procedure can change the relative weight of different NMEs in the fc^-factorization approach, 
that is out of our present study. 

Now we turn to the the -0(2 S) polarization issue, which is the most interesting part of 
our study. In general, the spin density matrix of a vector particle decaying into a lepton pair 
depends on three angular parameters Xg, X$ and A g^ which can be measured experimentally. 
The double differential angular distribution of the decay leptons can be written as [51] : 

—--—— ~ 1 + Xe cos 2 9* + Xs sin 2 6* cos 20* + Xgs sin 29* cos 0*, (19) 

d cos 9* dcj)* 9 ^ r 

where 9* and 0* are the polar and azimuthal angles of the decay lepton measured in the 
charmonium rest frame. The case of (Xg, X^, Xg$) = (0, 0, 0) corresponds to unpolarized state, 
while (Xg, Xf/,, Xg<j,) = (1, 0, 0) and (Xg, Xq, Xg$) = (—1, 0, 0) refer to fully transverse and longi¬ 
tudinal polarizations. The CMS [20] and LHCb [21] Collaborations have measured all these 
parameters as functions of 0(2 S) transverse momentum in two complementary frames: the 
Collins-Soper and hclicity ones. In addition, the CMS Collaboration provided measurements 
in the perpendicular hclicity frame. In the Collins-Soper frame the polarization axis z bisects 
the two beam directions whereas the polarization axis in the hclicity frame coincides with the 
direction of 0(2 S) momentum in the laboratory frame. In the perpendicular hclicity frame 
the z axis is orthogonal to that in the Collins-Soper frame and lies in the plane spanned by 
the two beam momenta. Additionally, the frame-independent polarization parameter [511152] 
A* = (Ag + 3A^)/(I — A 0 ) was investigated. Below we estimate the polarization parameters 
Xg, X 0 , A 00 and A* for the CMS and LHCb conditions. Our calculation generally follows the 
experimental procedure. We collect the simulated events in the kinematical region defined by 
the CMS and LHCb experiments, generate the decay lepton angular distributions according 
to the production and decay matrix elements, and then apply a three-parametric fit based 
on (19). 

In Figs. 4 — 8 we confront our predictions for polarization parameters Xg, X^, Xg^ and 
A* with the latest CMS [20] and LHCb [21] data. We find slight transverse polarization 
(A 0 ~ 0.2) in the Collins-Soper frame and slight longitudinal polarization (Xg ~ —0.2) in the 
hclicity frame at low transverse momenta covered by the LHCb experiment. These results 




are practically independent on the 0(2 S) rapidity. At higher p t the polarizations of '0(25') 
mesons, calculated in the Collins-Soper frame, go from slight transverse (A g ~ 0.15) to 
almost zero values (Xg ~ 0.05) as the transverse momentum increases from p? ~ 10 GeV to 
50 GeV. In the helicity and perpendicular helicity frames 0(25) polarization changes from 
longitudinal (A# ~ —0.2) to slight longitudinal (Ag ~ —0.1). Here we arrive at the key point 
of our paper. Figs. 4 — 8 show that 0(2 S) production, calculated in the /^-factorization 
approach, tends to be unpolarized at high pr, in agreement with the CMS data m- Indeed, 
as a strict consequence of the initial gluon off-shellness, a large fraction of 0(2 S) mesons 
with zero helicity is produced in the partonic subprocesses, including both CS and CO 
contributions. Moreover, the fraction of such events increases when p t grows up (see, for 
example, [27ll34] i. The only exception refers to 1 S^ channel, which is produced unpolarized 
due to its spinless nature. It is a remarkable property of the kr -factorization scheme that the 
gluon fragmentation to 3 S[ 8 ^ states produces nearly unpolarized 0(2 S) mesons (in contrast 
with conventional collinear NRQCD where the mesons carry strong transverse polarization). 
Thus, we can conclude that the problem of 0(2 S) spin alignment can be solved if the initial 
gluon off-shellness is taken into account. 

A comparison of our predictions with the LHC data [201121] shows that the latter seem to 
support the trend observed in the /^-factorization formalism. However, while our predictions 
for A^ and A g^ parameters agree with the data, the description of A g and A* is still rather 
qualitative than quantitative, due to the huge experimental uncertainties. 

Finally, we would like note that there are significant theoretical uncertainties connected 
with the choice of the renormalization and/or factorization scales, the inclusion of NLO 
subprocesses and exact definition of NMEs. The detailed study of these uncertainties is out 
of our present paper. 


4 Conclusions 


We have considered prompt 0(2 S) production and polarization in pp collisions at the LHC 
energy yfs = 7 TeV in the framework of /^-factorization approach. We have used the LO non- 
relativistic QCD formalism including both color-singlet and color-octet contributions. Using 
the TMD gluon densities in a proton derived from the CCFM equation and from the Kimber- 
Martin-Ryskin prescription, we extracted the color-octet NMEs (O^ 


ig( 8 ) 


and {OX [ 3 P 0 (8) ) for 0(25) mesons from fits to transverse momentum distributions provided 
by the latest CMS measurements. Using the fitted NMEs, we have successfully described 
the data presented by the ATLAS, CMS and LHCb Collaborations. We estimated the 
polarization parameters A g, A^ and A g^ and demonstrated that taking into account the off- 
shellness of the initial gluons in the color-octet contributions leads to unpolarized 0(25') 
production at high transverse momenta, that is in qualitative agreement with the LHC 
data. 
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Figure 1: The double differential cross sections of prompt ij}(2S) meson production at the 
LHC. Left panel: the dashed, dash-dotted, dotted and short dash-dotted curves correspond 
to the color-singlet 3 S'| 1 ' ) and color-octet 1 Sq 8 \ 3 S[ 8 \ 3 Pj 8 ' 1 contributions calculated with the 
KMR gluon density. The solid curve represent the sum of CS and CO terms. Right panel: 
the solid, dashed and dash-dotted curves correspond to the predictions obtained with the AO, 
JH and KMR gluon distributions, respectively. The experimental data are from ATLAS [16]. 
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Figure 2: The transverse momentum distribution of prompt ip(2S) meson production in pp 
collisions at the LHC. Notation of all curves is the same as in Fig. 1. The experimental data 
are from CMS [17] . 




Figure 3: The transverse momentum distribution of prompt ip (2S) meson production in pp 
collisions at the LHC. Notation of all curves is the same as in Fig. 1. The experimental data 
are from LHCb [IB] . 
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Figure 4: Polarization parameters Xg, A^, A g^ and A* of prompt mesons calculated as 

a function of if) (2S) transverse momentum in the Collins-Soper frame. The solid, dashed and 
dash-dotted histograms correspond to the predictions obtained at \y\ < 0.6, 0.6 < \y\ < 1.2 
and 1.2 < \y\ < 1.5. The KMR gluon distribution is used. The experimental data are from 
CMS [20]. ^ 
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Figure 5: Polarization parameters Xg, A^, Xg# and A* of prompt mesons calculated 

as a function of '0(2*5) transverse momentum in the helicity frame. The solid, dashed and 
dash-dotted histograms correspond to the predictions obtained at \y\ < 0.6, 0.6 < \y\ < 1.2 
and 1.2 < \y\ < 1.5. The KMR gluon distribution is used. The experimental data are from 
CMS [20]. ^ 
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Figure 6: Polarization parameters Xg, A^, A g^ and A* of prompt mesons calculated as 

a function of -0 (2S') transverse momentum in the perpendicular helicity frame. The solid, 
dashed and dash-dotted histograms correspond to the predictions obtained at \y\ < 0.6, 
0.6 < | y\ < 1.2 and 1.2 < \y\ < 1.5. The KMR gluon distribution is used. The experimental 
data are from CMS [20]. 
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Figure 7: Polarization parameters Xg, X^, A^ and A* of prompt '0(25') mesons calculated as a 
function of 0(25) transverse momentum in the Collins-Soper frame. The solid, dashed, dash- 
dotted, dotted and short dash-dotted histograms correspond to the predictions obtained at 
2 < y < 2.5, 2.5 < y < 3, 3 < y < 3.5, 3.5 < y < 4 and 4 < y < 4.5. The KMR gluon 
distribution is used. The experimental data are from LHCb [21]. 
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Figure 8: Polarization parameters Xg, X^, Xg^ and A* of prompt 0(2 S ) mesons calculated as 
a function of 0(2 S) transverse momentum in the helicity frame. The solid, dashed, dash- 
dotted, dotted and short dash-dotted histograms correspond to the predictions obtained at 
2 < y < 2.5, 2.5 < y < 3, 3 < y < 3.5, 3.5 < y < 4 and 4 < y < 4.5. The KMR gluon 
distribution is used. The experimental data are from LHCb [21] . 
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